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ABSTRACT 

The DNA-binding domains (DBDs) of class I steroid 
receptors— androgen, glucocorticoid, progesterone 
and mineralocorticoid receptors— recognize a 
similar c/s-element, an inverted repeat of 5 -AGAA 
CA-3 with a 3-nt spacer. However, these receptors 
regulate transcription programs that are largely 
receptor-specific. To address the role of the DBD 
in and of itself in ensuring specificity of androgen 
receptor (AR) binding to chromatin in vivo, we 
used SPARKI knock-in mice whose AR DBD has 
the second zinc finger replaced by that of the gluco- 
corticoid receptor. Comparison of AR-binding 
events in epididymides and prostates of wild-type 
(wt) and SPARKI mice revealed that AR achieves se- 
lective chromatin binding through a less stringent 
sequence requirement for the 3 -hexamer. In par- 
ticular, a T at position 12 in the second hexamer is 
dispensable for wt AR but mandatory for SPARKI 
AR binding, and only a G at position 11 is highly 
conserved among wt AR-preferred response 
elements. Genome-wide AR-binding events agree 
with the respective transcriptome profiles, in that 
attenuated AR binding in SPARKI mouse epididymis 
correlates with blunted androgen response in vivo. 
Collectively, AR-selective actions in vivo rely on 
relaxed rather than increased stringency of cis- 
elements on chromatin. These elements are, in 
turn, poorly recognized by other class I steroid 
receptors. 



INTRODUCTION 

Androgens mediate the effects of male sex steroids 
through the androgen receptor (AR) that regulates expres- 
sion of androgen-dependent genes by binding to cis- 
elements in the regulatory regions of target genes. The 
DNA-binding domain (DBD) of class I steroid recep- 
tors — AR, glucocorticoid receptor (GR), progesterone 
receptor (PR) and mineralocorticoid receptor (MR) — is 
highly conserved. They all recognize a similar palindromic 
response element that comprises an inverted repeat of the 
5'-AGAACA-3' hexamer with a 3-nt spacer, usually 
termed as a canonical androgen/glucocorticoid response 
element (ARE/GRE) (1). The first hexamer is highly 
conserved and represents the high-affinity site for 
receptor interaction because mutations in this core, as 
opposed to the second one, have a strong impact on 
receptor binding and reporter gene responsiveness, as 
shown by in vitro assays (2). In transient transfection ex- 
periments, AR, GR, PR and MR are all capable of 
binding to the canonical ARE/GRE and inducing expres- 
sion of reporter genes driven by the canonical sequence 
(3). Moreover, genome-wide analyses using chromatin 
immunoprecipitation (ChlP) coupled with high- 
throughput sequencing (ChlP-seq) experiments have 
shown that AR and GR bind to shared loci on native 
chromatin (4-6). One-half of the AR cistrome overlaps 
with that of GR in LNCaP-lF5 prostate cancer cells, 
and ligand-occupied GR is able to modulate the AR 
pathway (5). These results present a biological enigma: 
how are specific responses to cognate ligands of class I 
steroid receptors ensured in cells that express multiple re- 
ceptors at the same time. 
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One potential explanation for specific responses to class 
I steroid receptor ligands is subtle differences in the 
cM-elements that receptor DBDs recognize on chromatin. 
Selective AREs resembling direct repeats of the 5'-AGAA 
CA-3' hexamer are not recognized by GR or MR in vitro 
(3,7,8). A chimeric AR, in which the second zinc finger of 
the DBD is replaced by that of GR, possesses significantly 
attenuated affinity for selective AREs in transient trans- 
fection experiments and in vitro binding assays, whereas 
transactivation by the chimeric receptor via the classical 
ARE/GREs is similar to or better than that by wild-type 
(wt) AR (9). A transgenic mouse line with the same AR to 
GR switch in the second zinc finger of the DBD, dubbed 
SPARKI (specificity-affecting AR knock-in) mice, 
has smaller reproductive organs and reduced fertility 
compared with wt littermates, although androgen-depend- 
ent anabolic parameters, such as body weight and muscle 
strength, are unaffected (9). SPARKI males show reduced 
sperm count and impaired sperm maturation, suggesting 
defective function of the epididymis, a highly androgen- 
responsive reproductive organ (10). 

SPARKI mice provide an excellent in vivo model system 
to investigate the importance of DBD structure in and of 
itself in setting apart genome-wide androgen-specific re- 
sponses from those potentially regulated by glucocortic- 
oids as well. Here, we compared AR-binding events and 
transcription programs in epididymides and prostates of 
wt and SPARKI male mice to delineate potential genome- 
wide rules behind androgen-selective responses. Our 
results show that, counter-intuitively, AR selectivity is 
achieved by relaxed rather than tightened c/s-element 
stringency at chromatin binding sites which, in turn, at- 
tenuates binding of other class I steroid receptors to these 
cw-elements. 

MATERIALS AND METHODS 

Mice and hormone treatments 

All the experiments were carried out using 12- week-old adult 
male mice of C57BL/6 strain. SPARKI knock-in mice have 
been described previously (9). Wt littermates were used as 
controls. For ChIP assays, caput epididymides and ventral 
prostates were collected from intact wt and SPARKI male 
mice. For mRNA isolation, caput epididymides were har- 
vested from castrated wt and SPARKI mice. Testosterone 
treatment was brought about by sc Silastic implants (Silclear 
Tubing; Degania Silicone, Jordan Valley, Israel) for 7 days, 
corresponding to a daily dose of 23 (ig. Control group 
received empty implants. All experimental work involving 
animals was conducted with approval of the Katholieke 
Universiteit Leuven ethical committee. 

Chromatin immiinoprecipitation 

ChIP assays were performed as described previously (4) with 
modification for tissues (11). Ventral prostates from four 
and caput epididymides from six wt mice and eight 
SPARKI mice were pooled to yield one 500 -ul chromatin 
sample. Frozen tissues were pulverized and cross-linked with 
1% formaldehyde for 20min at room temperature. 
Formaldehyde was quenched by adding 0.125 M glycine 



for 5min at room temperature, followed by two washes 
with ice-cold phosphate-buffered saline. The cross-linked 
tissues were homogenized in hexylene glycol buffer (1 M 
hexylene glycol, 0.1 mM MgCi2, 5mM ethylene glycol 
tetraacetic acid, 1 mM PIPES, 2 mM dithiothreitol, 1 x pro- 
teinase inhibitor cocktail; Roche Diagnostics, Indianapolis, 
IN, USA) using Ultra Turrax (Ika, Staufen, Germany), 
filtered through nylon net (20 um), and nuclei isolated by 
low-speed centrifugation. The crude nuclear pellet was re- 
suspended in radioimmunoprecipitation assay buffer and 
sonicated using a micro-tip sonicator (Misonix Inc., 
Farmingdale, NY, USA) to yield chromatin fragments of 
100-500 bp in size. Anti-AR (12) and normal rabbit IgG 
(sc-2027, Santa Cruz, Dallas, TX, USA) antibodies were 
conjugated to Dynabeads (Invitrogen, Carlsbad, CA, 
USA) and incubated with sonicated chromatin overnight 
at 4°C, followed by five washes in LiCl buffer and incubation 
at 65°C overnight for reverse crosslinking. 
Immunoprecipitated DNA was purified and used for quan- 
titative polymerase chain reaction (qPCR) analysis. In 
ChlP-qPCR, immunoprecipitated and input DNA were 
amplified using SYBR Green Mastermix (Roche) and 
specific primers (Supplementary Table SI), and the results 
are shown as percentage of input values. 

ChlP-seq library preparation 

ChIP samples were processed according to Illumina's 
library preparation protocol by pooling three immunopre- 
cipitates for each library (4). In brief, DNA samples were 
blunt-ended, A-tailed and ligated to sequencing adapters 
followed by size selection (size range 150-300 bp) on 
agarose gel. Excised fragments were purified using 
Qiaquick Gel Extraction kit (Qiagen GmbH, Hilden, 
Germany) followed by PCR amplification for 20 cycles. 
The purified DNA was sequenced using GAII (Illumina 
Inc., San Diego, CA, USA) at Biomedicum Functional 
Genomics core facility. ChlP-seq reads were filtered ac- 
cording to Illumina's instructions using the Illumina 
chastity filter, and reads were aligned to mouse genome 
(mm9) using Bowtie without any mismatches. ChlP-seq 
experiments were carried out in biological duplicates. 

Bioinformatics analysis 

ChlP-seq peak calling and differential peak calling were 
performed using MACS2 (13,14) and HOMER (15) algo- 
rithms. The overlap analysis, tag density maps and 
binding site correlation plots were performed using 
Cistrome (16). The de novo motif analysis was performed 
using MEME Suite (17) and HOMER motif discovery 
algorithm (15). Data visualization was carried out using 
Integrative Genomics Viewer (18). Gene expression 
profiles for differentially expressed genes in whole epididy- 
mis of intact wt and SPARKI mice (10) were compared 
with the differentially androgen-regulated genes in caput 
epididymis (11). 

Electrophoretic mobility shift assay 

AR and GR DBD fragments of ~11 kDa in size, corres- 
ponding to residues 533-637 of rat AR DBD and residues 
432-533 of rat GR DBD, were expressed as glutathione 
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S-transferase fusion proteins in the Escherichia coli 
BL21 strain (8). The GST tag was removed by 
thrombin cleavage. Double-stranded oligonucleotides 
containing ARE sequences (Supplementary Table SI) 
were labeled with [oc- 32 P]dCTP by a filling-in with the 
Klenow fragment of DNA polymerase. Slp-MUT and 
Slp-HRE sequences have been described (3). Purified 
AR DBD or GR DBD (lOOng protein) was incubated 
with radiolabeled probe (20000cpm) in lOmM HEPES 
(pH 7.9), 2.5 mM MgCl 2 , 0.05 mM ethylenediaminete- 
traacetic acid, 8% glycerol, 50 mM NaCl, 2.5ng/ul 
poly (dldC), ImM dithiothreitol and 0.05% Triton 
X-100 for 20min on ice. DBD-bound and unbound 
probes were separated by electrophoresis on a 4% 
polyacrylamide gel. Radioactive bands were visualized 
with STORM 840 Phosphorlmager (Molecular 
Dynamics). 

Transactivation assays 

HeLa cells (ATCC, Manassas, VA, USA) were cultured 
in Dulbecco's modified Eagle's medium containing 4.5 g/1 
glucose, 4mM L-glutamine, penicillin (100IU/ml)- 
streptomycin (lOOug/ml) (Sigma-Aldrich, St. Louis, 
MO, USA) and 10% fetal calf serum (Invitrogen). For 
transactivation experiments, cells were seeded onto 96- 
well plates (8000 cells/well) in Dulbecco's modified 
Eagle's medium supplemented with 5% charcoal- 
stripped serum (Sigma-Aldrich) and transfected using 
GeneJuice (Novagen, Merck KGaA, Darmstadt, 
Germany) with lOng of AR or GR expression plasmid 
(19), 100 ng of luciferase reporter plasmid containing 
four copies of the ARE of interest (Supplementary 
Table SI) and 5ng of (3-galactosidase expression 
plasmid (Stratagene, Santa Clara, CA, USA). Reporter 
plasmids were generated as previously described (3). Cells 
were lysed after a 24-h exposure to lOnM 
methyltrienolone (Perkin Elmer, Waltham, MA, USA) 
or lOnM dexamethasone (Sigma-Aldrich). Luciferase 
and p-galactosidase activities were measured, and trans- 
activation activity was calculated (20). 

RNA isolation and quantitative reverse transcription PCR 

Caput epididymides of wt and SPARKI mice (« = 5 in 
both cases) were harvested in RNALater (Qiagen) after 
a 7-day treatment with testosterone or vehicle. Total 
RNA was isolated using RNeasy Mini Kit (Qiagen). 
cDNA synthesis was carried out from 2 ug of total RNA 
with Transcriptor High Fidelity cDNA synthesis kit using 
random hexamers (Roche). Quantitative reverse 
transcription PCR (qRT-PCR) was performed using 
SYBR Green mastermix (Roche) and normalized to 18 S 
rRNA levels. Student's f-test (P < 0.001) was used to cal- 
culate the statistical significance of differences in gene ex- 
pression in wt and SPARKI mice between testosterone- 
and vehicle-treated groups and between testosterone- 
treated wt and SPARKI groups. The primer sequences 
are provided in Supplementary Table SI. 



RESULTS 

Genome-wide AR-binding events in wt and SPARKI 
epididymis 

ChlP-seq analyses of in vivo AR-binding sites (ARBs) on 
epididymal chromatin were performed in biological dupli- 
cate samples. Examples of AR ChlP-seq tracks of wt and 
SPARKI mice are shown in Figure 1A at two genomic loci 
that are androgen-responsive in murine epididymis (11); 
Fkbp5 (FK506 binding protein 5) encoding a chaperone 
protein and Acox3 (acyl-coenzyme A oxidase 3) involved 
in fatty acid metabolism in peroxisomes. AR occupancy is 
indistinguishable at the Fkbp5 regulatory region between 
wt and SPARKI mice, which agrees with the result that 
Fkbp5 mRNA accumulation is induced by androgen to the 
same extent in epididymides of these mice (10). Sequenced 
reads were aligned to mouse genome (mm9) using bowtie, 
and peak calling for individual replicates was carried out 
using MACS2 (13,14). Importantly, the biological repli- 
cates show excellent concordance, as judged by individual 
ChlP-seq tracks (Figure 1A) and by coefficient values 
from correlation analyses of the raw sequence data 
(Figure IB). To increase the depth of analysis, we 
concatenated independent biological replicates and re-per- 
formed peak calling using two different algorithms, 
MACS2 and HOMER (15). The peaks called from the 
two algorithms were used to compute the final overlapping 
peaks (with at least 1 nt overlap), resulting in AR 
cistromes comprising 10009 and 6446 ARBs in wt and 
SPARKI epididymides, respectively (Figure 1C, 
Supplementary Data sets SI and S2). 

SPARKI AR fails to recognize wt AR-preferred enhancers 

Two-thirds of the SPARKI AR cistrome overlap with that 
of wt mice (Figure 2A). The greater number of ARBs that 
are unique to wt epididymis is explainable by the ability of 
wt AR to recognize selective AREs. Because the SPARKI 
AR is anticipated to recognize the canonical ARE/GRE 
only, it is intriguing that there are also in vivo AR-binding 
events that are favored by the SPARKI receptor 
(Figure 2A). As ChlP-seq provides a quantitative assess- 
ment of receptor binding across the genome, the tag 
density of AR binding to a given site could be used as a 
surrogate measure for the affinity of the receptor to this 
site. The tag density maps of AR-binding events and 
average tag profiles in wt and SPARKI mice indicate 
that the two receptors are loaded to the same extent 
onto sites shared by wt AR and SPARKI AR, most 
likely reflecting indistinguishable binding affinities of the 
receptors for these sites (Figure 2B and C). Overall, wt and 
SPARKI ARs bind to the shared sites with ~2-fold higher 
affinity than to non-shared sites (Figure 2C). Interestingly, 
shared sites have also been reported to be the high-affinity 
binding loci for the estrogen receptor, when chromatin 
binding events in multiple estrogen target tissues or cell 
lines were compared with those with strict tissue- or cell- 
type specificity (21). 

The 'uniqueness' of wt and SPARKI AR-binding events 
is relative rather than absolute, in that there are significant 
quantitative differences in the recruitment of the two ARs 
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Figure 1. In vivo AR binding in epididymis of wt and SPARKI mice. (A) AR-binding events mapped to Fkbp5 and Acox3 loci in caput epididymis of 
wt and SPARKI male mice. ChlP-seq tracks from two independent biological replicates and non-specific rabbit IgG control are shown. (B) Genome- 
scale wiggle file correlation plot and correlation coefficient values for AR-binding events of two independent biological replicates in epididymides of 
wt and SPARKI mice. (C) The number of AR-binding peaks derived from concatenated biological replicates and using MACS2 and HOMER peak 
calling algorithms. (FDR, false discovery rate). 



to these sites. Therefore, we call the two ARB groups as wt 
AR-preferred and SPARKI AR-preferred sites (Figure 2B 
and C). The tag density maps constructed for the wt AR- 
preferred sites with a > 4-fold and a 2^1-fold change over 
SPARKI ARBs indicate that wt AR-binding events with 
high fold changes are poorly recognized by SPARKI AR 
(Figure 2B). The reverse is true for the AR-binding events 
preferred by SPARKI AR; these sites are recognized less 
well by wt AR (Figure 2B). The average tag profiles — a 
more quantitative assessment of the differences — reveal 
that the ARBs preferred by wt AR have tag numbers 
three times higher than those of SPARKI AR at the 
same sites, and the binding events preferred by SPARKI 
AR have twice the number of tags over those of wt AR at 
the same sites (Figure 2C). Collectively, these data indicate 
that there are classes of c«-elements that wt and SPARKI 
ARs recognize with dissimilar binding affinities. 

Examples of ChlP-seq tracks are shown in Figure 2D 
with sequences at the summit of a shared AR-binding 
event ( 5' -acAGA AC AgaaTGTTCTct- 3') and a SPARKI 
AR-preferred site (5'-agAGAACAtggTGTTCCct-3')- 
Both possess a canonical ARE/GRE, a palindrome 



of 5'-AGAACA-3' with a 3-nt spacer. The example of a 
cw-element preferred by wt AR is a non-canonical 
ARE (Figure 2D) that bears resemblance to an imperfect 
direct repeat sequence (5'-gaGGAACAgggTGAACC 
aa-3'). 

WT AR selectivity in vivo depends on relaxed cis-element 
stringency 

We hypothesized that genome-wide AR-binding events in 
epididymides of wt and SPARKI mice will help in estab- 
lishment of defined rules for the cw-elements that AR rec- 
ognizes on chromatin. To this end, we performed de novo 
motif search analyses using MEME Suite (17) for the three 
ARB categories: (i) binding sites shared by the two recep- 
tors or (ii) preferred by wt AR or (iii) SPARKI AR — on 
sequences comprising the top 1000 sites spanning ± 15 bp 
around the summit of the receptor binding locus. As 
expected, an inverted palindromic repeat with a 3-nt gap 
resembling the canonical ARE/GRE is highly enriched 
among the shared ARBs (E = 3.7e-639). Interestingly, a 
similar cw-element is also highly enriched among the 
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Figure 2. Features of the AR-binding events in epididymides of wt and SPARKI mice. (A) Area-proportional Venn diagrams of AR cistromes in 
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ARBs preferred by SPARKI AR (E = 3.2e-437). In both 
instances, there is a high degree of sequence conservation. 
The best conserved nucleotides are a G at position 2 and a 
C at position 5, followed by an A at positions 4 and 6, and 
their complementary partners T, G, T and C at positions 
10, 11, 12 and 14, respectively (Figure 3A). 

De novo motif search analyses of the binding sites 
preferred by wt AR uncovered a different cw-element 
(E = l.le-403), in which the 5' hexamer is almost identical 
with that described above but separated by a 3-nt gap 
from a hexamer with weak sequence conservation 
(Figure 3A). This latter hexamer has a highly conserved 
G at position 11, along with weak conservation for T and 
C at positions 10 and 14, respectively. Importantly, the 
highly conserved T at position 12 of the canonical ARE/ 
GRE present at the shared and SPARKI AR-preferred 
sites is dispensable among the c«-elements for wt AR 
binding. Thus, relaxed stringency in the c;,v-element 
sequence — dismissal of the conserved T at position 12 in 
particular — ensures selective wt AR binding to chromatin. 
This result agrees with a previous mutational analysis 
showing that substitution of an A with a T at this 



position in a selective ARE enhances GR binding and 
glucocorticoid response (2). 

Wt AR-preferred and SPARKI AR-preferred binding 
events were next cross-compared directly by HOMER. 
This analysis revealed that 31% of the wt AR-preferred 
peaks have >4-fold enrichment over those to SPARKI 
AR (1721 ARBs, Supplementary Data set S3), and 21% 
of the SPARKI AR-preferred peaks have >4-fold enrich- 
ment over those of wt AR (414 ARBs, Supplementary Data 
set S4) (Figure 3B). Importantly, <1% of the shared AR- 
binding peaks show >4-fold differential enrichment, 
implying an equal binding affinity for wt and SPARKI 
AR at the shared sites (Figure 3B). To characterize 
further the nature of the cw-elements among the wt AR- 
preferred ARBs, de novo motif search analysis was per- 
formed on the top 500 ARBs with >4-fold difference 
between wt and SPARKI AR binding peaks (Figure 3B, 
uppermost and lowest lines). This analysis yielded binding 
motifs for wt and SPARKI AR that are similar to those 
obtained by direct de novo motif searches (cf. Figure 3A 
and C). Of note, different analysis parameters, such as dif- 
ferent windows around the peak summits, the use of 
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Figure 3. Sequence enrichment analysis and validation of the AR-binding events. (A) Cu-elements identified by de novo motif search for the ARBs 
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shuffled instead of top-enriched sites or different algorithms 
(HOMER or MEME Suite) did not change the results of 
de novo analysis. Direct ChlP-qPCR validation of 13 
randomly chosen loci corresponding to AR-binding 
events in the three categories confirm that there are signifi- 
cant quantitative differences in the loading of AR onto sites 
preferred either by wt AR or SPARKI AR (Figure 3D). 

In vitro validation of the wt AR-preferred, SPARKI 
AR-preferred and shared cis-elements 

Electrophoretic mobility shift assay (EMSA) with AR 
DBD and GR DBD was used to examine the specificity 
of a few ARE sequences in the three categories under 
in vitro conditions. Classical and selective AREs, 
characterized in previous work (3), were used as 
controls. The second zinc finger of SPARKI AR is iden- 
tical with that of GR, and SPARKI AR-preferred and 
shared sequences were bound by both AR DBD and 
GR DBD (Figure 4A and B). The AREs classified by 
our ChlP-seq analysis as wt AR-preferred sequences 
were not recognized by GR DBD at all or with low 
affinity (Figure 4C). Thus, the results from EMSAs 
using naked DNA and AR DBD or GR DBD agree 
with those of ChlP-seq experiments carried out in a 
genuine chromatin environment and with full-length 
receptors. 

To validate the ARB categories further, three AREs 
from each group were cloned in four copies upstream of 



the luciferase reporter gene and exposed in transactivation 
assays in HeLa cells to methyltrienolone and dexametha- 
sone for AR and GR, respectively. AR was able to 
transactivate via all nine AREs, confirming their func- 
tional enhancer nature (Figure 4D). Importantly, GR 
could transactivate only through shared and SPARKI- 
preferred ARBs comprising canonical AREs but not 
through wt AR-preferred ARE sequences (Figure 4E). 

Differential AR-binding events are linked to differentially 
expressed genes in wt and SPARKI epididymis 

By examining the gene expression data from whole 
epididymides of intact wt and SPARKI male mice (10), 
we found 5 1 1 transcripts with higher expression in wt than 
SPARKI mice ( = wt-preferred genes) and 843 transcripts 
with higher expression in SPARKI than wt mice ( = 
SPARKI-preferred genes) (fold change > 1.5, P<0.05). 
Even though intact male mice have physiological 
androgen levels, these differentially expressed transcripts 
may not all be androgen-regulated. To address this issue, 
we compared the Kerkhofs et al. data set (10) with 
androgen-regulated genes in caput epididymides of 
castrated male mice treated for 3 days with testosterone 
or vehicle (11) and identified 219 androgen-regulated tran- 
scripts that are differentially expressed in intact SPARKI 
versus wt mice (110 androgen upregulated and 109 
androgen downregulated transcripts, Supplementary 
Data set S5). Unsupervised hierarchical clustering of the 
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Figure 4. Validation of the three ARB categories by EMSAs and transactivation assays. Radiolabeled double-stranded oligonucleotides containing 
the AREs were tested for binding by purified AR DBD or GR DBD. Arrows refer to free probe and arrowheads to ARE DBD complexes. Slp-MUT 
and Slp-HRE are classical and selective ARE controls, respectively. (A) AREs from three SPARKI AR-preferred binding sites. (B) AREs from three 
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or dexamethasone (Dex) (E) for 24 h. Results are normalized to P-galactosidase activities and presented as mean + SEM values of three biological 
replicates performed in triplicate. 



219 transcripts revealed that most of the wt-pref erred 
genes are upregulated by androgen, whereas most of the 
SPARKl-preferred genes are downregulated by androgen 
exposure (Figure 5A). Importantly, the wt-pref erred 
androgen-regulated genes such as Rhbg, Rhbdd3, Ehd3, 
Isynal, Acat2, Ccno, Ramp3, Lcnl2 and Cyp4fl5 show 
blunted upregulation by testosterone treatment in 
SPARKI epididymis compared with wt, as judged by 
mRNA accumulation (Figure 5B). Moreover, loss or at- 
tenuation of AR binding at regulatory sites of the wt- 
pref erred genes Rhbg, Rhbdd3, Ehd3, Isynal, Acat2 and 
Ccno is observed in SPARKI epididymis compared with 
wt (Figure 5C). Likewise, qRT-PCR assays confirm the 
differential expression of SPARKl-preferred genes in 
SPARKI epididymis compared with the wt epididymis 
as seen in the case of Arrdc2, Pdk2ll and CalmB 
(Figure 6A). There is greater loading of SPARKI AR 
than wt AR onto the Pkd2ll locus, but in some instances, 
such as Arrdc2, there is no clear difference between wt AR 
and SPARKI AR- binding events (Figure 6B). Examples 
of androgen-regulated genes (Fkbp5, Acox3, Lipg, Creld2 
and Kcnkl) that are stably expressed in both wt and 
SPARKI epididymides are shown in Figure 6C. In 



agreement with transcript accumulation, AR loading is 
equal onto the respective regulatory regions in wt and 
SPARKI epididymis, as exemplified by Fkbp5 and 
Acox3 (Figure 1A) and Lipg and Creld2 (Figure 6D). 

Relaxed AR c/s-element stringency is not limited to the 
epididymis 

To examine whether the finding on in vivo AR selectivity 
also applies to another androgen target tissue, we per- 
formed AR ChlP-seq on ventral prostates of SPARKI 
mice and their wt littermates. Peak calling with MACS2 
and HOMER identified 6693 and 5552 high-confidence 
AR-binding events in prostates of wt and SPARKI 
mice, respectively, 3370 of which are shared by the two 
receptors (Figure 7A, Supplementary Data sets S6 and 
S7). SPARKI AR has low affinity to wt AR-preferred 
sites and wt AR to SPARKI AR-preferred sites, whereas 
both receptors bind to shared sites with high affinity 
(Figure 7B). Of the wt AR-preferred ARBs, 1079 
sites have >4-fold enrichment over SPARKI AR 
(Supplementary Data set S8), and 609 SPARKI 
AR-preferred sites have >4-fold enrichment over wt AR 
(Supplementary Data set S9). Importantly and similar to 
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Figure 5. Comparison of AR-binding events and transcriptome profiles in wt and SPARKI epididymis. (A) Unsupervised hierarchical clustering of 
transcripts that are regulated by testosterone in caput epididymis (T versus vehicle) and expressed differentially in epididymides of intact SPARKI 
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white bar (vehicle). (C) ChlP-seq tracks of six representative examples showing differences in AR loading onto genomic loci adjacent to AR- 
regulated genes in wt and SPARKI epididymides. 



epididymis, the cw-element identified by de novo motif 
search among the wt AR-preferred ARBs in prostate 
shows weak conservation of the second hexamer 
sequence, highlighting the relaxed ra-element stringency, 
with the G at position 1 1 being an important determinant 
for selective AR binding in vivo (Figure 7C). In a manner 
identical with that in epididymis, SPARKI AR-preferred 
and shared ARBs are enriched in prostate for the canon- 
ical ARE/GRE (Figure 7C). Examples of AR ChlP-seq 
tracks at three androgen-responsive loci in murine 
prostate, Pgap2, Creld2 and Slc25a30 (11) depict wt and 
SPARKI AR loading onto enhancers of the three ARB 
categories — wt and SPARKI AR-preferred or shared — 
and show the respective ARE sequences at peak 
summits (Figure 7D). 



DISCUSSION 

Previous work has shown that there are two types of cis- 
elements that are recognized by the AR: (i) classical AREs 
corresponding to the canonical ARE/GRE and recognized 



by other class I steroid receptors as well, and (ii) selective 
AREs recognized primarily by AR and showing features 
that resemble direct repeats of the 5'-AGAACA-3' hexamer 
(3,10,22). By using genome-wide in vivo AR ChlP-seq on 
epididymal and prostatic chromatin of wt and SPARKI 
mice, we show in this work that the SPARKI AR fails to 
bind or binds with blunted affinity to sites corresponding to 
selective AREs. A significant proportion — 55 and 50% of 
the epididymal and prostatic AR cistromes, respectively — 
comprised wt AR-preferred ARBs affirming the biological 
importance of selective AREs in ensuring specific AR re- 
sponses. Surprisingly, the genome-wide analyses revealed 
that the selectivity of AR binding over that of other class 
I steroid receptors relies on relaxed rather than more strin- 
gent requirement for the response element sequence. 
Attenuated binding of SPARKI AR to selective AREs 
was, in turn, accompanied by blunted response to andro- 
gens of the corresponding genes in SPARKI mice. 

The 5' hexameric core is almost identical between canon- 
ical ARE/GREs and selective AREs, and this hexamer is 
suggested to be the high-affinity DNA sequence for the 
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binding of the first DBD (22,23). The contact of the first 
DBD to the high-affinity site is assumed to lead to a con- 
formational change in the receptor protein, and possibly also 
in the DNA element, which subsequently facilitate binding 
of the second DBD to the less conserved low-affinity site 
(22,23). Previous work has indicated that the 3' hexamer 
of selective AREs resembles, in many instances, an imperfect 
direct repeat of the 5' hexameric sequence 5'-AGAACA-3' 
[summarized in (22)]. Our current results show why this 
would be the case, in that the Gs at positions 2 and 1 1 are 
highly conserved in the selective AREs, permitting 3' 
hexameric sequences that can be interpreted as being imper- 
fect direct repeats. More importantly, however, selective 
AREs exhibit no requirement for a T at position 12. 
Genome-wide ChlP-seq experiments on several cell lines or 
tissues have shown that m-elements for GR and PR contain 
a highly conserved T at this position (4—6,24—27). Thus, AR 
selectivity is achieved by loosening the structural stringency 
in the 3' hexameric sequence, and consequently, other class I 
steroid receptors fail to bind to this element. Crystal struc- 
ture studies have indicated that AR DBD has substitutions 
that create higher affinity for the receptor dimer interface 
than that of GR DBD, and subsequently, the AR could bind 



to a more diverse set of ris-elements with higher affinity and 
specificity than the GR (22,23). Thus, the high-affinity 
protein-protein binding at dimer interface rather than strin- 
gent cw-element structure is important for AR selectivity 
over the other class I type receptors. Other features in the 
receptor structure could possibly also contribute to the dif- 
ferent binding affinities of AR and GR as well, as mutations 
changing the amino acids in the dimerization interface did 
not affect in vitro DNA-binding affinities or transactivation 
capabilities of the two receptors (28). 

We examined in this work the role of DBD in and of 
itself in assuring the selectivity of AR signaling by 
comparing wt and chimeric receptors. However, specific 
transcriptional outcome in vivo by steroid receptors is 
influenced not only by the ra-element and its interaction 
with the DBD, but also by a number of chromatin modi- 
fications and divergent collaborating factors (29). Pioneer 
factors are proteins capable of binding to compact chro- 
matin and priming it for steroid receptor binding. In 
prostate cancer cells and in normal mouse prostate, AR 
binding is primed by FoxAl, and a distinct m-element, 
composed of an ARE half-site and a FoxAl element, is 
highly enriched among the AR-binding sequences (4,11). 



Nucleic Acids Research, 2014, Vol. 42, No. 7 4239 



Wild-type AR 

6,693 
sites 



B 



in 

m 

on 

< 




SPARKI AR 

5,552 
sites 



15- 




\ ■ wt 






\ • SPARKI 


10- 






5- 








-400 


0 400 




2 
a. 
a 

E> 15- 



a> 



-400 


0 400 




/ \ ■ wt 
\ ■ SPARKI 


-400 


1 — 
0 


400 



Distance from center (bp) 



Wild-type over SPARKI AR 

E = 3.0e-404 

2.0- 



Wild-type AR-preferred 



i i i i i i 



i i i i i i i i 

2 4 6 8 10 12 14 



Shared 

E = 2.8e-334 



r 11 
igG|_ 

r i! 

wt AR 
SPARKI AR [ 1 



Chr7:109,372kb 
_l 



109,384 kb 
I 



Pgap2 



AAACACT T CAAC 



Shared 




[21 

igG 



Chr1 5:88,652 kb 
_l 



wt AR | 
SPARKI AR [ 



I I I I I I 

2 4 6 8 10 12 14 



1 



Creld2 



ACAC T TTCT 



SPARKI over wild-type AR 

E= 1.5e-174 



SPARKI AR-preferred 



2.0 -, 

l %ii-LlT 



i i i i i i 

2 4 6 



IgG I 

noi 

wt ar| 

moi 

SPARKI AR 



Chr14:76,180 kb 
I 



8 10 12 14 



Slc25a30 



A AACA AAT TCCT 



Figure 7. AR cistromes in wt and SPARKI prostates. (A) Area-proportional Venn diagrams of AR cistromes in prostates of wt and SPARKI mice. 
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Pioneer factors have been reported to guide lineage-, cell- 
type- and tissue-specific AR binding (5,11,30), and the 
interplay between steroid receptors and pioneer factors 
is likely to contribute to the specificity of steroid 
receptor binding in addition to the DNA-DBD inter- 
action. However, we did not find markedly dissimilar 
compilations of m-elements for other transcription 
factors adjacent to the three ARE categories that would 
explain their divergent receptor binding properties 
(Supplementary Data set S10). 

The determinants of preferential in vivo binding of 
SPARKI AR to a subset of ARBs remain elusive. The 
DNA sequence in itself may act as an allosteric ligand 
for nuclear receptors, thereby affecting receptor conform- 
ation and activity (31). However, it is unlikely that this 
phenomenon is the main reason for the presence of 
SPARKI AR-preferred ARBs in vivo, as the SPARKI 
AR sites contain the same ARE sequence as the shared 
sites. Moreover, the cw-elements of SPARKI AR- 
preferred ARBs in isolation, i.e., in EMSA experiments 
and reporter gene assays (cf. Figure 4), lost their selectivity 
and responded equally well to wt AR and GR, supporting 
the importance of local chromatin environment. 



Differences in local chromatin structure, such as histone 
and DNA modifications, could explain the occurrence of 
these sites (32). Finally, in addition to cw-elements, cross- 
talk of steroid receptors with tissue-specific pioneer factors 
and/or coregulators is likely to influence receptor recruit- 
ment to selected chromatin sites and determine the final 
transcriptional outcome. 

ACCESSION NUMBERS 

The ChlP-seq data have been deposited in the Gene 
Expression Omnibus database with accession number 
GSE51106, and the gene expression data with accession 
number GSE47181. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 
ACKNOWLEDGEMENTS 

Saija Kotola and Rita Bollen are thanked for excellent 
technical assistance. 



4240 Nucleic Acids Research, 2014, Vol. 42, No. 7 



FUNDING 

Academy of Finland, Sigrid Juselius Foundation, Finnish 
Cancer Foundations, Biocentrum Helsinki, Helsinki 
University Central Hospital, Helsinki Biomedical 
Graduate Program and research grant of the KU 
Leuven nr OT/ 11/081. Holder of a doctoral fellowship 
of the Fund for Scientific Research Flanders (FWO) (to 
V.D.). Funding for open access charge: Sigrid Juselius 
Foundation and Biocentrum Helsinki. 

Conflict of interest statement. None declared. 



REFERENCES 

1. Mangelsdorf,D.J., Thummel,C, Beato.M., Herrlich,P., Schtttz,G, 
Umesono,K., Blumberg.B., Kastner,P., Mark,M., Chambon,P. 

et al. (1995) The nuclear receptor superfamily: the second decade. 
Cell, 83, 835-839. 

2. Verrijdt.G., Schoenmakers.E., Haelens.A., Peeters,B., 
Verhoeven.G., Rombauts.W. and Claessens.F. (2000) Change of 
specificity mutations in androgen-selective enhancers. Evidence for 
a role of differential DNA binding by the androgen receptor. 

J. Biol. Chem., 275, 12298-12305. 

3. Denayer.S., Helsen.C, Thorrez.L., Haelens.A. and Claessens.F. 
(2010) The rules of DNA recognition by the androgen receptor. 
Mol. Endocrinol, 24, 898-913. 

4. Sahu.B., Laakso.M., Ovaska.K., Mirtti.T., LundinJ., Rannikko.A., 
Sankila.A., TurunenJ.P., Lundin,M., KonstiJ. et al. (2011) Dual 
role of FoxAl in androgen receptor binding to chromatin, androgen 
signalling and prostate cancer. EMBO J., 30, 3962-3976. 

5. Sahu.B., Laakso.M., Pihlajamaa,P., Ovaska,K., SinielnikovJ., 
Hautaniemi,S. and Janne.O.A. (2013) FoxAl specifies unique 
androgen and glucocorticoid receptor binding events in prostate 
cancer cells. Cancer Res., 73, 1570-1580. 

6. John,S., Sabo,P.J., Thurman,R.E., Sung,M.H., Biddie,S.C, 
Johnson,T.A., Hager,G.L. and StamatoyannopoulosJ.A. (2011) 
Chromatin accessibility pre-determines glucocorticoid receptor 
binding patterns. Nat. Genet., 43, 264-268. 

7. Claessens.F., Alen,P., Devos,A., Peeters.B., Verhoeven,G. and 
Rombauts,W. (1996) The androgen-specific probasin response 
element 2 interacts differentially with androgen and glucocorticoid 
receptors. /. Biol. Chem., 271, 19013-19016. 

8. Schoenmakers.E., Alen.P., Verrijdt,G, Peeters,B., Verhoeven,G., 
Rombauts.W. and Claessens,F. (1999) Differential DNA binding 
by the androgen and glucocorticoid receptors involves the second 
Zn-finger and a C-terminal extension of the DNA-binding 
domains. Biochem. J.. 341, 515-521. 

9. Schauwaers.K., De Gendt.K., Saunders.P.T., Atanassova,N., 
Haelens.A., Callewaert.L., Moehren.U., SwinnenJ.V., 
Verhoeven.G, Verrijdt,G. et al. (2007) Loss of androgen receptor 
binding to selective androgen response elements causes a 
reproductive phenotype in a knockin mouse model. Proc. Natl 
Acad. Sci. USA, 104, 4961^1966. 

10. Kerkhofs.S., Dubois,V., De Gendt,K., Helsen,C, Clinckemalie,L., 
SpansX., Schuit,F., Boonen.S., Vanderschueren.D., Saunders.P.T. 
et al. (2012) A role for selective androgen response elements in 
the development of the epididymis and the androgen control of 
the 5alpha reductase II gene. EASES J., 26, 4360-4372. 

1 1. Pihlajamaa,P., Sahu,B., Lyly,L., Aittomaki.V., Hautaniemi,S. and 
Janne.O.A. (2014) Tissue-specific pioneer factors associate with 
androgen receptor cistromes and transcriptional programs. 
EMBO J., 33, 312-326. 

12. Kang.Z., Janne.O.A. and PalvimoJ.J. (2004) Coregulator 
recruitment and histone modifications in transcriptional regulation 
by the androgen receptor. Mol. Endocrinol., 18, 2633-2648. 

13. Feng, J., Liu,T., Qin,B., Zhang,Y. and Liu,X.S. (2012) Identifying 
ChlP-seq enrichment using MACS. Nat. Protoc, 7, 1728-1740. 

14. Zhang,Y., Liu,T., Meyer,C.A., EeckhouteJ., Johnson,D.S., 
Bernstein,B.E., Nusbaum.C, Myers, R.M., Brown,M., Li,W. et al. 



(2008) Model-based analysis of ChlP-Seq (MACS). Genome Biol, 
9, R137. 

15. Heinz.S., Benner.C, Spann,N, Bertolino,E., Lin.Y.C, Laslo,P., 
ChengJ.X., Murre.C, Singh.H. and Glass.C.K. (2010) Simple 
combinations of lineage-determining transcription factors prime 
cis-regulatory elements required for macrophage and B cell 
identities. Mol. Cell, 38, 576-589. 

16. Liu.T., OrtizJ.A., Taing,L., Meyer.C.A., Lee,B., Zhang,Y., 
Shin,H., Wong.S.S., Ma,J., Lei,Y. et al. (2011) Cistrome: an 
integrative platform for transcriptional regulation studies. 
Genome Biol., 12, R83. 

17. Bailey,T.L., Boden,M., Buske.F.A., Frith,M., Grant.C.E., 
Clementi.L., Ren,J., Li,W.W. and Noble.W.S. (2009) MEME 
SUITE: tools for motif discovery and searching. Nucleic Acids 
Res., 37, W202-W208. 

18. RobinsonJ.T., Thorvaldsd6ttir,H., Winckler.W., Guttman,M., 
Lander,E.S., Getz,G. and Mesirov,J.P. (2011) Integrative 
genomics viewer. Nat. Biotechnol., 29, 24-26. 

19. Scheller.A., Hughes,E., Golden,K.L. and Robins.D.M. (1998) 
Multiple receptor domains interact to permit, or restrict, 
androgen-specific gene activation. /. Biol. Chem., 273, 
24216-24222. 

20. Verrijdt.G, Schauwaers.K., Haelens,A., Rombauts.W. and 
Claessens.F. (2002) Functional interplay between two response 
elements with distinct binding characteristics dictates androgen 
specificity of the mouse sex-limited protein enhancer. /. Biol. 
Chem., 277, 35191-35201. 

21. GertzJ., Savic.D., Varley,K.E., Partridge,E.C, Safi,A., Jain,P., 
Cooper,G.M., Reddy,T.E., Crawford,G.E. and Myers.R.M. (2013) 
Distinct properties of cell-type-specific and shared transcription 
factor binding sites. Mol. Cell, 52, 25-36. 

22. Helsen,C, Kerkhofs,S., Clinckemalie,L., Spans,L., Laurent,M., 
Boonen,S., Vanderschueren.D. and Claessens.F. (2012) Structural 
basis for nuclear hormone receptor DNA binding. Mol. Cell 
Endocrinol, 348, 411-417. 

23. Shaffer.P.L., Jivan,A., Dollins.D.E., Claessens,F. and 
Gewirth,D.T. (2004) Structural basis of androgen receptor binding 
to selective androgen response elements. Proc. Natl Acad. Sci. 
USA, 101, 4758^1763. 

24. Ballare,C, Zaurin,R., Vicent.G.P. and Beato,M. (2013) More help 
than hindrance: nucleosomes aid transcriptional regulation. 
Nucleus, 4, 189-194. 

25. Rubel.C.A., Lanz,R.B., Kommagani,R., Franco, H.L., Lydon.J. 
and DeMayo.F.J. (2012) Genome-wide profiling of progesterone 
receptor binding in the mouse uterus. Mol. Endocrinol., 26, 
1428-1442. 

26. Lain,A.R., Creighton.C.J. and Conneely,O.M. (2013) Progesterone 
receptor targetome underlying mammary gland branching 
morphogenesis. Mol. Endocrinol., 27, 1743-1761. 

27. Grontved,L., John.S., Baek.S., Liu.Y., BuckleyJ.R., Vinson.C, 
Aguilera.G. and Hager.G.L. (2013) C/EBP maintains chromatin 
accessibilty in liver and facilitates glucocorticoid receptor 
recruitment to steroid response elements. EMBO J., 32, 
1568-1583. 

28. Verrijdt.G, Tanner.T., Moehren.U., CallewaertX., Haelens,A. 
and Claessens,F. (2006) The androgen receptor DNA-binding 
domain determines androgen selectivity of transcriptional 
response. Biochem. Soc. Trans., 34, 1089-1094. 

29. Zaret.K.S. and CarrollJ.S. (2011) Pioneer transcription factors: 
establishing competence for gene expression. Genes Dev., 25, 
2227-2241. 

30. Lupien.M., EeckhouteJ., Meyer,C.A., Wang,Q., Zhang,Y., Li,W., 
CarrollJ.S., Liu,X.S. and Brown,M. (2008) FoxAl translates 
epigenetic signatures into enhancer-driven lineage-specific 
transcription. Cell, 132, 958-970. 

31. Meijsing.S.H., PufalhM.A., So.A.Y., Bates,D.L., Chen,L. and 
Yamamoto,K.R. (2009) DNA binding site sequence directs 
glucocorticoid receptor structure and activity. Science, 324, 
407^410. 

32. Wiench,M., Miranda.T.B. and Hager,G.L. (2011) Control of 
nuclear receptor function by local chromatin structure. FASEB J., 
278, 2211-2230. 



